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Abstract 

A nonperturbative renormalization of the t^** model is considered. First we inte- 
grate out only a single pair of conjugated modes with wave vectors ±q. Then we are 
looking for the RG equation which would describe the transformation of the Hamil- 
tonian under the integration over a shell A — dA < fc < A, where dA 0. We show 
that the known Wegner-Houghton equation is consistent with the assumption of a 
simple superposition of the integration results for ±q. The renormalized action can 
be expanded in powers of the v?'' coupling constant u in the high temperature phase 
at M — J> 0. We compare the expansion coefficients with those exactly calculated by the 
diagrammatic perturbative method, and find some inconsistency. It causes a question 
in which sense the Wegner-Houghton equation is really exact. 

1 Introduction 

The renormalization group (RG) approach, perhaps, is the most extensively used one in 
numerous studies of critical phenomena [H [2] . Particularly, the perturbative RG approach 
to the or Ginzburg-Landau model is widely known O HI [U [6]. However, the pertur- 
bative approach suffers from some problems [?]■ Therefore it is interesting to look for a 
nonperturbative approach. Historicahy, nonperturbative RG equations have been devel- 
oped in parallel to the perturbative ones. These are so cahed exact RG equations (ERGE). 
The method of deriving such RG equations is close in spirit to the famous Wilson's ap- 
proach, where the basic idea is to integrate out the short-wave fluctuations corresponding 
to the wave vectors within A/s < g < A with the upper (or ultraviolet) cutoff parameter 
A and the renormalization scale s > 1. The oldest nonperturbative equation of this kind, 
originally presented by Wegner and Houghton [8], uses the sharp momentum cutoff. Later, 
a similar equation with smooth momentum cutoff has been proposed by Polchinski [9]. The 
RG equations of this class are reviewed in [TO] . 

According to the known classification |101 [TT] , there is another class of nonperturbative 
RG equations proposed by [12] and reviewed in [11]. Some relevant discussion can be found 
in |10] . as well. Such equations describe the variation of an average effective action rfc[0] 
depending on the running cutoff scale k. Here i;^(x) = ((/'(x)) is the averaged order- 
parameter field (for simplicity, we refer to the case of scalar field) . According to [12] , the 
averaging is performed over volume ~ k~'^ such that the fiuctuation degrees of freedom 
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with momenta q > k are effectively integrated out. In fact, the averaging over volume 
~ k~'^ is the usual block-spin-averaging procedure of the real-space renormalization. 
At the same time, the fluctuations with q < k are suppressed by a smooth infrared 
cutoff. As one can judge from [TTj, the existence of a deterministic relation between the 
configuration of external source {J(x)} and that of the averaged order parameter {(/>(x)} is 
(implicitly) assumed in the nonperturbative derivation of the RG flow equation. Namely, 
it is stated (see the text between (2.28) and (2.29) in [TT]) that 5J{'x)/5(j){y) is the inverse 
of (5<?!»(x)/5J(y), which has certain meaning as a matrix identity. To make this point 
clearer, let us consider a toy example J = Acj), where J = ( J(xi), J(x2), . . . , J(xn)) and 
(/) = ((ji(xi), (/'(X2), . . . , 0(xn)) are A'^-component vectors and A is a matrix of size N x N. 
In this case 9J(xj)/9(/>(xj) is the element Aij of matrix A, whereas d(p{xi)/dJ{xj) is the 
element (^~^).. of the inverse matrix A~^. In the continuum limit A^ — )• 00, this toy 
example corresponds to a linear dependence between { J(x)} and {(/){x)}. The calculation 
of derivative always implies the linearisation around some point, so that the matrix identity 
used in pA] (as a continuum limit in the above example) has a general meaning. However, 
it makes sense only if there exists a deterministic relation between the configurations of 
0(x) and J(x) or, in a mathematical notation, if there exist mappings / : { J(x)} — )■ {(/)(x)} 
and : {(/>(x)} — {J(x)}. On the other hand, according to the block-averaging, the 
values of should be understood as the block-averages. These, of course, are not 
uniquely determined by the external sources, but are fluctuating quantities. So, we are 
quite sceptical about the exactness of such an approach of averaged effective action. 

The integration over fluctuation degrees with momenta q > k does not alter the be- 
havior of the infrared modes, directly related to the critical exponents. From this point of 
view, the approach based on the equations of Wegner-Houghton and Polchinski type seems 
to be more natural. These are widely believed to be the exact RG equations, although, in 
view of our currently presented results, it turns out to be questionable in which sense they 
are really exact. In any case, the nonperturbative RG equations cannot be solved exactly, 
therefore a suitable truncation is used. The convergence of several truncation schemes and 
of the derivative expansion has been widely studied in [HI El [l5l [l6l [IT] . Here [IT] refers 
to the specific approach of [12]. A review about all the methods of approximate solution 
can be found in [T8]. 

Another problem is to test and verify the nonperturbative RG equations, comparing 
the results with the known exact and rigorous solutions, as well as with the results of 
the perturbation theory. In |15j , the derivative expansion of the RG /3-function has been 
considered, showing the agreement up to the second order between the perturbative results 
and those obtained from the Legendre flow equation, which also belongs to the same class 
of RG equations as the Wegner-Houghton and Polchinski equations. It has been stated 
in [13] that the critical exponent v, extracted from the Wegner-Houghton equation in the 
local potential approximation, agrees with the e-expansion up to the 0{e) order, as well 
as with the 1/n (1/A^ in the notations of |13] ) expansion in the leading order. However, 
looking carefully on the results of |13j . one should make clear that "the leading order of 
the 1/n expansion" in this case is no more than the zeroth order, whereas the expansion 
coefficient at 1/n is inconsistent with that proposed by the perturbative RG calculation 
at any fixed dimension d except only d = 4. The inconsistency could be understood from 
the point of view that the Wegner-Houghton equation has been solved approximately. 
Therefore it would be interesting to verify whether the problem is eliminated beyond the 
local potential approximation. One should also take into account that the perturbative 
RG theory is not rigorous and, therefore, we think that a possible inconsistency still would 



not prove that something is really wrong with the nonperturbative RG equation. In any 
case, it is a remarkable fact that correct RG eigenvalue spectrum and critical exponents 
are obtained in the local potential approximation at n — )• oo from the Wegner-Houghton 
equation [13], as well as from similar RG equations [19], in agreement with the known 
exact and rigorous results for the spherical model. It shows that some solutions, being 
not exact, nevertheless can lead to exact critical exponents. From this point of view, it 
seems also possible that some kind of approximations, made in the derivation of an RG 
equation, are not harmful for the critical exponents. 

We propose a simple test of the Wegner-Houghton equation: to verify the expansion of 
the renormalized action S of the ip^ model in powers of the coupling constant u at u — )• 
in the high-temperature phase. Such a test is rigorous, in the sense that the natural 
domain of validity of the perturbation theory is considered. We think that it would be 
quite natural to start with such a relatively simple and straightforward test before passing 
to more complicated ones, considered in [131 [El IS]- We have made this simplest test 
in our paper and have found that the Wegner-Houghton equation fails to give all correct 
expansion coefficients. We have also proposed another derivation of the Wegner-Houghton 
equation (Sees. [H [3]). It is helpful to clarify the origin of the mentioned inconsistency. It 
is also less obscure from the point of view that the used assumptions and approximations 
are clearly stated. As regards the derivation in [8], at least one essential step is obscure 
and apparently contains an implicit approximation which, in very essence, is analogous to 
that pointed out in our derivation. We will discuss this point in Sec. [3l 

2 An elementary step of renormalization 

To derive a nonperturbative RG equation for the model, we should start with some 
elementary steps, as explained in this section. 

Consider the action S[ip] which depends on the configuration of the order parameter 
field (/j(x) depending on coordinate x. By definition, it is related to the Hamiltonian H of 
the model via S = H/T, where T is the temperature measured in energy units. In general, 
(p{x) is an n-component vector with components <^j(x) given in the Fourier representation 
as ipj{x) = Sfc<A where V = L'^ is the volume of the system, d is the spatial 

dimensionality, and A is the upper cutoff of the wave vectors. We consider the action of 
the Ginzburg-Landau form. For simplicity, we include only the ip"^ and ip'^ terms. The 
action of such 93^ model is given by 

S[^] ='^®(X)^j,k^j-k + uV''^ X] '/'i,klV'i,k293/,k3V'/ -ki-k2-k3 , (1) 
i,k j,«,ki,k2,k3 

where 0(k) is some function of wave vector k, e. g., 0(k) = tq + ck'^ like in theories of 
critical phenomena [H E] El [7] . In the sums we set 99; = for A; > A. 

The renormalization group (RG) transformation implies the integration over (/jj ^ for 
some set of wave vectors with A' < /c < A, i. e., the Kadanoff's transformation, followed 
by certain rescaling procedure [1]. The action under the Kadanoff's transformation is 
changed from S[ip] to Stra.['p] according to the equation 




Alternatively, one often writes —Stra[^] + AL instead of — S'tralv'] to separate the constant 
part of the action AL'^. This, however, is merely a redefinition of Stra; and for our purposes 



it is suitable to use Q. Note that (^j^k = V'j k + ^V'j'k ^ complex number and = 
holds (since ^Pj{x) is always real), so that the integration over ipj^^ means in fact the 
integration over real and imaginary parts of (/7j^k for each pair of conjugated wave vectors 
k and — k. 

The Kadanoff's transformation ([2]) can be split in a sequence of elementary steps 
S[ip] — )• S'tralv'] of the repeated integration given by 



-S'tra[</3] 



(3) 



for each j and q G where Q is the subset of independent wave vectors (itq represent 
one independent mode) within A' < g < A. Thus, in the first elementary step of renor- 
malization we have to insert the original action ([1]) into ([3]) and perform the integration 
for one chosen j and q € 1^. In an exact treatment we must take into account that the 
action is already changed in the following elementary steps. 

For A' > A/3, we can use the following exact decomposition of ([T]) 



5[(/7] =Ao+ ^l(/7j,q + ^IV^i-q + ^2¥'j,q'/'j -q + ^aV'j, q + ^2V'j -q + ^4V'j, q^^j-q 



where 
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Here the sums are marked by a prime to indicate that terms containing 9?j,±q are omitted. 
This is simply a splitting of ([T]) into parts with all possible powers of 'Pj,±ci- The condition 
A' > A/3, as well as the existence of the upper cutoff for the wave vectors, ensures that 
terms of the third power are absent in Besides, the derivation is performed formally 
considering all (/?;^k as independent variables. 

Taking into account , as well as the fact that Ai = A'^ + iA'( and B2 = B2 + iB!^ are 
complex numbers, the transformed action after the first elementary renormalization step 
reads 



Straif] = Aq 




exp [-2 q - ^'lV"q)] 



-A2 [ + f 



^4 U'j.q +^[ 



(10) 



Considering only the field configurations which are relevant in the thermodynamic 
limit V oo, Eq. (jlOp can be simplified, omitting the terms with B2 and A4. Really, 
using the coordinate representation ipn^ = V^^l"^ J e~^^^dx, we can write 



Bo = 2uV-^ 



^(1 + 25^,) I (/,f(x)e^^-dx 



) ■ (11) 



The quantity f ipf{'x.)e'^'^^dx. is an average of ^fi'x) over the volume with oscillating 
weight factor e**'^. This quantity vanishes for relevant configurations in the thermody- 
namic limit: due to the oscillations, positive and negative contributions are similar in 
magnitude and cancel at F — > 00. Since (| P) = (| P) is the Fourier transform 

of the two-point correlation function, it is bounded at y — )• 00 and, hence, V^^-q also 
is bounded for relevant configurations giving nonvanishing contribution to the statistical 
averages (•) in the thermodynamic limit. Consequently, for these configurations, A2 is 
a quantity of order 0{1), whereas V~^ip'j_^ and B2 vanish at y — )• 00. Note, however, 
that the term with A^ = O {V~^^ in (jlOp cannot be neglected unless A2 is positive. One 
can judge that the latter condition is satisfied for the relevant field configurations due to 
existence of the thermodynamic limit for the RG fiow. 

Omitting the terms with B2 and A4, the integrals in (jlOp can be easily calculated. It 
yields 

StrM = S'M+ASt[^], (12) 

where S'[ip] = Aq is the original action, where only the itq modes of the j-th field com- 
ponent are omitted, whereas AS'f^'^[(/7] represents the elementary variation of the action 
given by 

A5.1bl=ln(ii)-^. (13) 

According to the arguments provided above, this equation is exact for the relevant field 
configurations with > in the thermodynamic limit. 

The contributions to ^ and ([7]) provided by modes with wave vectors k, obeying two 
relations A — dA < k < A and k 7^ ±q, are irrelevant in the thermodynamic limit at 
dA 0. It can be verified by the method of analysis introduced in Sec. [2j Hence, Eq. (jl3p 
can be written as _ _ 

AS^iM = In ( 4^ ) - + 5Sty] , (14) 




where 



~ dS A-dA 

^1 = -P^— =4^^"^ V^j,ki</5i,k2y'« -q-ki-ka , (15) 

A2 = Pj: ^ = G(q)+G(-q)+4nF-i V (l + 2,5zj) I 99z,k |2 , (16) 



i,k 



and bSf^^ip\ is a vanishingly small correction in the considered limit. Here the operators 
P set to zero all (^j^k within the shell A — dA < A; < A (i. e., the derivatives are evaluated 
at zero 'Pj^ for k within the shell), and the upper border A — dA for sums implies that we 
set ffiy = for k > A — dA. The above replacements are meaningful, since they allow to 
obtain easily the Wegner-Houghton equation, as discussed in the following section. 



3 Superposition hypothesis and the 
Wegner— Houghton equation 

Intuitively, it could seem very reasonable that the result of integration over Fourier modes 
within the shell A — dA < /c < A at dA — )• can be represented as a superposition of 
elementary contributions given by (jl4p . neglecting the irrelevant corrections (5S't/jj [93] . We 
will call this idea the superposition hypothesis. 

We remind, however, that strictly exact treatment requires a sequential integration of 
exp{—S[ip]) over a set of fj^q- The renormalized action changes after each such integration, 
and these changes influence the following steps. A problem is to estimate the discrepancy 
between the results of two methods: (1) the exact integration and (2) the superposition 
approximation. Since it is necessary to perform infinitely many integration steps in the 
thermodynamic limit, the problem is nontrivial and the superposition hypothesis cannot 
be rigorously justified. 

Nevertheless, the summation of elementary contributions in accordance with the su- 
perposition hypothesis leads to the known Wegner-Houghton equation [8j. In this case 
the variation of the action due to the integration over shell reads 



It is exactly consistent with Eq. (2.13) in [5]. The factor 1/2 appears, since only half of 
the wave vectors represent independent modes. Here we have indicated that the quantities 



figuration [ip]. If Ai and A2 are represented by the derivatives of S[(p] (see ([T5]) and (fT6]) ). 
then the equation is written exactly as in [8]. 



To avoid possible confusion, one has to make clear that the operators P in (|15p and (jl6p 
influence the result, as discussed further on. It means that the equation where these 
operators are simply omitted, referred in the review paper [10] as the Wegner-Houghton 
equation, is not really the Wegner-Houghton equation. 

The derivation in [8j is somewhat different. Instead of performing only one elementary 
step of integration first, the expansion of Hamiltonian in terms of all shell variables is made 
there. The basic method of [8] is to show that, in the thermodynamic limit at dA 0, the 
expansion consists of terms containing no more than two derivatives with respect to the 
field components. Moreover, it is assumed implicitly that only the diagonal terms with 
k' = — k are important finally, when performing the summation over the wave vectors 
k, k'. It leads to Eq. (2.12) in [8]. The omitting of nondiagonal terms is equivalent to the 
superposition assumption we discussed already. Indeed, in this and only in this case the 
integration over the shell variables can be performed independently, as if the superposition 
principle were hold. Hence, essentially the same approximation is used in [8] as in our 
derivation, although it is not stated explicitly. 

Our derivation refers to the 93^ model, whereas in the form with derivatives the equation 
may have a more general validity, as supposed in [8]. Indeed, ()14p remains correct for a 
generalized model provided that higher than second order derivatives of S[ip] vanish for 
relevant field configurations in the thermodynamic limit. It, in fact, has been assumed 
and shown in [8]. Based on similar arguments we have used already, the latter assumption 
can be justified for certain class of models, for which the action is represented by a linear 
combination of ip'^-kind terms with wave-vector dependent weights and vanishing sum of 




(17) 



Ai and A2 depend on the current j and q. They depend also on the considered field con- 



the wave vectors Ylh=i = related to the (f factors. In this case we have 

^l(j,q) = <^i,-q75 , (18) 

A2(j,q) = (19) 

for the relevant configurations at y — ?> oo and dA — > 0. The second term in (jlSh appears 
because the derivative dS/dipj^q contains relevant terms with v^j,-q, which have to be 
removed. The influence of the operators P is seen from (I15p and (llSp . 

Here we do not include the second, i. e., the rescaling step of the RG transformation. 
It, however, can be easily calculated for any given action, as described, e. g., in It is 
not relevant four our further considerations. 



4 The weak coupling limit 

Here we consider the weak coupling limit u — >■ of the model with B(k) = + ck^ at 
a given positive tq, i. e., in the high temperature phase. In this case A5tra[v] can be 
expanded in powers of u. It is the natural domain of validity of the perturbation theory, 
and the expansion coefficients can be calculated exactly by the known methods applying 
the Feynman diagram technique and the Wick's theorem (UEKTl]. On the other hand, the 
expansion can be performed in (jl7p . Our aim is to compare the results of both methods 
to check the correctness of (fT7|) . since the latter equation is based on assumptions. 

Let us denote by AStralv'] the variation of S[(p] omitting the constant (independent of 
the field configuration) part. Then the expansion in powers of u reads 

AStrM = ASi [^] u + [as':^^ [^] + A5f [ip] + AS':^^ [^]) u^ + O (n^) , (20) 

where the expansion coefficient at is split in three parts ASg"^ /S.S2^ [p], and IS.S^'^ [^p] 
corresponding to the ip^, ip'^, and tp^ contributions, respectively. The contribution of order 
u is related to the diagram O , whereas the three second-order contributions — to the 
diagrams — O — i XZX , and ^ — ^ . The diagram technique represents the 
expansion of —S[ip] in terms of connected Feynman diagrams, where the coupled lines are 
associated with the Gaussian averages. In particular, the Fourier transformed two-point 
correlation function in the Gaussian approximation Go(k) = {^j^k^j^-k)o = l/[2G(k)] 
appears due to the integration over (p'j^^ and y^Jj^. It is represented as the coupling of 
lines, in such a way that each line related to the wave vector k and vector-component j is 
coupled with another line having the wave vector — k and the same component j. Thus, if 
we integrate over 99^ ^ within A — dA < /c < A in ([2]) , then it corresponds to the coupling of 
lines in the same range of wave vectors in the diagram technique. According to the Wick's 
theorem, one has to sum over all possible couplings, which finally yields the summation 
(integration) over the wave vectors obeying the constraint A — dA < k < A for each of the 
coupled lines associated with the factors Go(k). In the n-component case, it is suitable 

to represent the ip'^ vertex as >• <C , where the same index j is associated with two 

solid lines connected to one node. The above diagrams are given by the sum of all possible 

couplings of the vertices >• <i , yielding the corresponding topological pictures when 

the dashed lines shrink to points. In this case factor n corresponds to each closed loop 
of solid lines, which comes from the summation over j. For a complete definition of the 



diagram technique, one has to mention that factors —uV~^ are related to the dashed 
Hnes, Go(k) - to the coupled solid lines, and the fields (pj^k " to the outer uncoupled solid 
lines. Besides, each diagram contains a combinatorial factor. For a diagram consisting 

of m vertices >• <C , it is the number of all possible couplings of (numbered) lines, 

divided by m!. 

At dA — )• 0, the diagrammatic calculation for the n~component case yields 

A-dA 

ASi^\ip] = -4V-^ '/'i,ki¥'i,k29'«,k3V'i-ki-k2-k3 (22) 

i,/,ki,k2,k3 

X [{n + 4)Q (ki + k2. A, dA) + iQ (ki + kg. A, dA)] 

A-dA 

AS^^'^Lp] = -8V~'^ ^ ¥'i,ki9'j,k2'/'j,k3V'i,k4'/'i,k5y'«,-ki-k2-k3-k4-k5 

j,j,J,ki,k2,k3,k4,k5 

xGo(ki + k2 + k3) ^(|ki+k2 + k3 \,A,dA) , (23) 

where = S{d)/{2^Y, S{d) = 2ti'^/'^ /T{d/ 2) is the area of unit sphere in d dimensions, 
0(A) is the value of 0(k) at A; = A, whereas J-{k, A, dA) is a cutoff function which has the 
value 1 within A — dA < k < A and zero otherwise. The quantity Q is given by 

Q(k,A,dA) = V-' Yl G'o(q)Go(k-q)^(|k-q|,A,dA) . (24) 

A-dA<q<A 



Below we will give some details of calculation of (|22p . which is the most important term 
in our further discussion. To obtain this result, we have dechipered the ^d^xC diagram 
as a sum of three diagrams of different topologies made of vertices > <C , i. e., 

>" — O ^ ) ^ C| , and I I , providing the same topological picture 

XZX when shrinking the dashed lines to points. Recall that any loop made of solid 

lines of >» <C gives a factor n, and one needs also to compute the combinatorial 

factors. For the above three diagrams, the resulting factors are 4n, 16, and 16, which 
enter the prefactors of Q in ()22p . To obtain the correct sign, we recall that the diagram 
expansion is for —S[ip]. The other diagrams are calculated in a similar way. 

The expansion of (I17p gives no contribution AS2°'\^], and we have skipped it in the 
diagrammatic calculation as an irrelevant term, which vanishes faster than oc dA at dA — )• 
in the thermodynamic limit V oo. The expansion of the logarithm term in (llTh yields 
AS'i[(^] exactly consistent with (j2ip . Similarly, A£'2'^^[(/9] is exactly consistent with (j23p . 

One has to remark that two propogators are involved in (j24p and, therefore, the volume 
of summation region with nonvanishing cut function shrinks as (dA)^ for a given nonzero 
wave vector k at dA — )• 0. However, there is a contribution linear in dA for k = 0. As a 
result, a contribution proportional to dA appears in ()22p . 

Note that the contributions ()2ip and (|23p come from diagrams with only one coupled 
line. The term ()22p is related to the diagram with two coupled lines. The expansion 
of p!7|) provides a different result for the corresponding part of AS'tra['/']: 

\ d—1 J A A—dA 

A^'^M = - Q2a) E + 4 + 45,7) 1 ? I ^^k2 P . (25) 



Note that (j25p comes from the In A2 term in (jl7p , and the calculation is particularly simple 
in this case, since the related sum in (jl6p is independent of q. Eq. (j25p is obtained if we 
set (5(k, A, dA) — )■ (^^,0 Q(0; A, dA) in (|22p (in this case only the diagonal terms j = I are 
relevant when summing up the contributions with Q (ki + ka, A, dA), as it can be shown by 
an analysis of relevant real-space configurations, since {ipj(x)(pi(x.)) = holds for j ^ I). 
It means that a subset of terms is missing in (j25p . as compared to (j22p . The following 
analysis will show that this discrepancy between (I22p and ()25p is important. 

It is interesting to mention that (|25p is obtained also by the diagrammatic perturbation 
method if we first integrate out only the mode with <Pj^±ci and then formally apply the 
superposition hypothesis, as in the derivation of the Wegner-Houghton equation. It shows 
that the discrepancy between (j25|) and (j22]) arises because in one case the superposition 
hypothesis is applied, whereas in the other case it is not used. 

The difference between (j22p and ()25p can be better seen in the coordinate representa- 
tion. In this case (I22p reads 



AS^^\ip]= - (An + IG) J J ip\xi)R'^{xi-X2)v\^2)dxidx2 (26) 

- 16^/ /"(/?j(xi)v?z(xi) i?^(xi - X2) 99j(x2)(/?i(x2) dxidx2 , 

where 

i?(x) = Go(q) T{q, A, dA)e'^^ (27) 

q 

is the Fourier transform of GqJ-, and (x) = '/'f (x). In three dimensions we have 

= (OT) r'"<'^^' <^«' 

for any given x at dA — t- and L 00, where L is the linear system size. The continuum 
approximation (I28p . however, is not correct for x ~ L and therefore, probably, should not 
be used for the evaluation of (j26p . 

The coordinate representation of (I25p is 



(6) i^dA'^-^dA 

M = 02(X^ 



(n + 4) / / ip\xi)V-^ip\x2)dj^idx2 



+4^ / / ¥''(xi)y-Vi(x2)c?xidx2 



(29) 



Eq. (j29p represents a relevant contribution at dA — t- 0, as it is proportional to dA. It is ob- 
viously not consistent with ()26p . In fact, the term (I29p represents a mean-field interaction, 
which is proportional to 1/V and independent of the distance, whereas ()26p corresponds 
to another non-local interaction given by i?^(xi — X2). Hence, the Wegner-Houghton 
equation (jlTp does not yield all correct expansion coefficients at n — )• 0. 



5 Discussion 



The results of our test, stated at the end of Sec. [H reveal some inconsistency between 
the Wegner-Houghton equation and the diagrammatic perturbation theory in the high 



temperature phase at u — > 0. Since this is the natural domain of vahdity of the per- 
turbation theory, there should be no doubts that it produces correct results here, which 
agree with ([2]). So, the results of our test point to some inconsistency between the Wegner- 
Houghton equation and ([2]) , which causes a question in which sense the Wegner-Houghton 
equation is really exact. The same can be asked about the equations of Polchinski type, 
since these (as it is believed) are generalizations of the Wegner-Houghton equation to 
the case of smooth momentum cutoff. There is no contradiction with the tests of consis- 
tency made in [13\ I15j. since our test is independent and quite different. According to our 
derivation of the Wegner-Houghton equation and the related discussion, it turns out that 
the reason of the inconsistency, likely, is the superposition approximation (defined at the 
beginning of Sec. [3]) used in our paper and analogous approximation implicitly used in [8]. 
Despite of this problem, the Wegner-Houghton equation is able to reproduce the exact 
RG eigenvalue spectrum and critical exponents of the spherical model at n — t- oo |13j . 
This fact can be interpreted in such a way that the superposition approximation (or an 
analogous approximation) is valid to derive such nonperturbative RG equations, which can 
produce correct (exact) critical exponents in some limit cases, at least. From a general 
point of view, it concerns the fundamental question about the relation between the form 
of RG equation and the universal quantities. It has been verified in several known studies 
that the universal quantities are invariant with respect to some kind of variations in the 
RG equation, like changes in the shape of the momentum cutoff function. This property 
is known as the reparametrisation invariance [10]. Probably, the universal quantities are 
invariant also with respect to such a variation of the Wegner-Houghton equation, which 
makes it exactly consistent with However, this is only a hypothesis. 

6 Conclusions 

1. The nonperturbative Wegner-Houghton RG equation has been rederived (Sees. [2] 
and[3l), discussing explicitly some assumptions which are used here. In particular, 
our derivation assumes the superposition of small contributions provided by elemen- 
tary integration steps over the short-wave fluctuation modes. We consider it as 
an approximation. As discussed in Sec. [3l the original derivation by Wegner and 
Houghton includes essentially the same approximation, although not stated explic- 
itly. 

2. According to our calculation in Sec. HI the Wegner-Houghton equation is not com- 
pletely consistent with the diagrammatic perturbation theory in the limit of small 
(p^ coupling constant u in the high temperature phase. This fact, together with 
some other important results known from literature, is discussed in Sec. O Apart 
from critical remarks, a hypothesis has been proposed that the equations of Wegner- 
Houghton type, perhaps, can give exact universal quantities. 
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